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ABSTRACT 



This paper presents and examines new near-infrared integral field observations of the three so-called 'embedded star clusters' lo- 
cated in the nuclear region of NGC1365. Adaptive-optics- corrected K-band data cubes were obtained with the ESO/VLT instrument 
SINFONI. The continuum in the K-band and emission lines such as He I, Bry, and several H2 lines were mapped at an achieved angular 
resolution of 0.2" over a field of 3x3 around each source. We find that the continuum emission of the sources is spatially resolved. 
This means that they are indeed cluster complexes confined to regions of about 50 pc extension. We performed robust measurements 
of the equivalent width of the CO absorption band at 2.3 pm and of Bry. For the main mid-infrared bright sources, the data only allow 
us to determine an upper limit to the equivalent width of the CO bands. Under the assumption of an instantaneously formed standard 
initial mass function Starburst99 model, the new measurements are found to be incompatible with previously published mid-infrared 
line ratios. We show that an upper mass limit of 25 to 30 Mq, lower than the typically assumed 100 Mq, allows one to simply remove 
this inconsistency. For such a model, the measurements are consistent with ages in the range of 5.5 Myr to 6.5 Myr, implying masses 
in the range from 3 to 10 x 10* Mq. We detect extended gas emission both in H 11 and Ht. We argue that the central cluster complexes 
are the sources of excitation for the whole nebulae, through ionisation and shock heating. We detect a blue wing on the Bry emission 
profile, suggesting the existence of gas outflows centred on the cluster complexes. We do not find any evidence for the presence of a 
lower mass cluster population, which would fill up a "traditional" power law cluster mass function. 
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1. Introduction 

The knowledge about young star cluster populations in galax- 
ies increased much since high angular resolution (~0.1") visi- 
ble images could be ob tained with the Hubble Space Telescope 
(see the reviews by Porteg ies Zwart et"aD 120101: iLarsenl 1201 ll 
and references therein). In the infrared, images approaching such 
a high angular resolution can be achieved only from large aper- 
ture, hence ground-based telescopes. The high angular resolu- 
tion infrared observations therefore suffer from the limitations 
inherent to the technique itself: only a few infrared transpar- 
ent windows exist, and, due to the intense atmospheric emis- 
sion, the sensitivity is relatively low. The few high angular res- 
olution mid-infrared (MIR) studies of actively star-forming re- 
gions in external galaxies have been able to resolve them in hot 
spots, which generally have no counterpart in the visible, but are 
bright in the radio. These types of sourc es have been name d "em- 
bedded cltisters'' (iGoriian et al.ll200ll: IpTa nte & Sauvag dl2002t 
iBeck et al."2002'; 'Marti'n-Hernandez et al. 2005; Galliano et al.' 
2005; Martin-Hernandez et al. 2006; Wold & Galliano 2006; 
Gandhi et al.l l20I ll) . They correspond to the hypothetical first 



phase of a star cluster, when it is still embedded in its birth ma- 
terial. It is believed that this embedded phase only lasts for a few 
million years: the intense stellar winds and supernova activity 
provokes the expulsion of the remaining gas and quick ly leads 
to the disruption of the cluster itself (see the review bv lBastianI 
I2OI iL and references therein). 



Three examples of massive star clusters in this an evolution- 
ary stage have been discovered in the star-fo rming nuclear regioii 
of the strongly barred spiral NGC1365 (Galliano et al. 20051 
GA05). The relatively small distance of NGC1365 (18.6Mpc, 
90 pc per arcsec), the fact t hat this galaxy h as been extensively 
observed and modelled (see ILindbladll 1 999 l for a review on the 
specific galaxy), and the fact that the galaxy hosts an active 
galactic nucleus (AGN) that can feed an adaptive optics (AO) 
system make these clusters prefeiTed targets for detailed studies. 
They are located in the central star-forming region, close to, or 
within, a well-defined dust lane. They were named M4, M5 and 
M6, 'M' standing for MIR (GA05). 

iGalliano et al.l (l2005h presented the discovery (with 
ESO/3.6m/TIMMI2) of the bright MIR counterparts to 
previously known centimetre sources dSandqvist et al. 1982 
ISaikia et al. 1994; San dqvist et al.lll995l: For bes & Norris 199^^ 
and in terpreted them as 'embedded star clusters'. . Galliano et al] 
(I2OO8I) . hereafter GA08, complemented these observations 
with near-infrared (NIR) data (K and L images and spectra 
with ESOA'LT/ISAAC), as well as new MIR data (narrow 
12.8 yum image and N-band spectra with ESO/VLT/VISIR). 
These sources are bright in the MIR continuum, nebular lines, 
cm continuum and ar e spatially associated to CO hot-spots 
(ISakamoto et all 12007'. hereafter SA07). In GA08, we deter- 
mined masses of about several times 10^ M© for the stellar 
component of each cluster, which put s them in the range of th e 
most massive known stellar clusters (iKissler-Patig et al]l2006l) . 
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The masses were estimated using the dereddened Bry fluxes, 
and a determination of the cluster ages. Because the hydrogen 
emission line fluxes dramatically decrease with time during 
the first ~10Myr of a cluster life, the mass estimated using 
the nebular line flux is very sensitive to the age determination. 
In GA08, we proposed an age determination for the clusters 
using the ionised gas emission lines in the N-band. Three 
lines of diff'erent ionisation potential ([Arlll] 27.6eV, [SIV] 
34.8eV and [Nell] 21.6eV) are located in the N-band. In the 
VISIR N-band spectra, only the lower ionisation line [Nell] 
could be detected. The stellar population synthesis model for 
young and instantaneously formed clusters predicts a drop of 
more than two orders of magnitude in the stellar continuum 
for radiation capable of ionising [Arlll] and [S IV], while this 
drop is much smaller (one order of magnitude) for the [Ne II] 
ionising radiation (See Fig. 15, 16 and 17 in GA08). From 
this argument, we estimated an age of about 7Myr for the 
three sources, from which we derived the above mentioned 
masses. Some other observational facts supported our age 
estimate for the clusters. First, similar-quality MIR spectra 
of a younger massive e mbedded cluster (NGC5253/C 2 with 
ESO/3.6m/TIMMI2 by iMartm-Hernandez et al.ll2005h indeed 
display intense [SIV]. They show that [SIV] is much brighter 
than [Nell]. Also, the centimetre spectral indexes of embedded 
clusters in NGC1365 reveal a significant share of non-thermal 
emission, possibly originating from a high supernovae rate. 
This points towards ages older than 3Myr or 4Myr (GA05). 
In contrast, you nger (~lMyr old) clusters h ave flat or even 
inverted spectra (iKobulnickv & Johnsonll 19991) . The long-slit K- 
band ISAAC spectra published in GA08 also displayed the CO 
absorption bands at 2.3 /im. These features are only expected 
to appear with the first red super giant stars, hence not before 
6Myr. Nevertheless, since the ISAAC spectra have relatively 
low angular resolution, this argument was not considered to 
be robust because we could not infirm the possibility that the 
absorption came from the background or nearby sources. A 
more reliable analysis is presented here, and interestingly, 
does not confirm the presence of the CO absorption bands that 
originates from the clusters themselves. 

Elmegreen et al.l (l2009h examined how these sources might 
have been formed. These authors showed that the extreme 
stellar masses derived for the clusters in GA08 are consis- 
tent with the properties of the galaxy: the large amounts 
of gas found in the central kpc regio n that was detected 
through HI jj orsater & van Moorsell[T995l) and CO observations 
(ISakamoto et al...2007|) are sufficient for the formation of a few 
10^ Mq clusters. lElmegreen et al.l (l2009t) propose a scenario to 
explain the building-up of the nuclear mass, where gas loops in 
along dust filaments to feed the bar dust lane and finally flows 
inwards towards the nuclear region. The interaction between the 
dust filaments and the dust lane is possibly a trigger for the for- 
mation of star clusters, which then would move inwards follow- 
ing the gas flow and feed the central bulge. 

In the present paper, we discuss a new dataset for these 
three targets, recorded with SINFONI, the adaptive optics near- 
infrared integral field instrument available at the VLT. The ob- 
servations were obtained with the adaptive optics loop closed on 
the nucleus of the galaxy, and hence have an angular resolution 
that is close to the diffraction limit of the 8m aperture telescope. 
Since the SINFONI data cubes contain very much information, 
we give priority to the detailed presentation of the data, includ- 
ing the publication of tables. We report about the collection and 
reduction of the data in Sec. |2] The data analysis including the 
different measurement procedures used are described in Sec. |3] 



In Sec. |4] the observational results are provided; they are dis- 
cussed in Sec.|5] and our conclusions are given in Sec.|6l 

2. Data collection and reduction 

The dataset was obtained with SINFONI at the ESOA'LT 
in the night from 2005/12/04 to 2005/12/Oil SINFONI is a 
NIR integral field spectrometer (SPIFFI) coupled to a visi- 
ble curvature AO system. The wavefront was analysed directly 
on the AGN of the galaxy (type 1.5, R=15.4). The foflow- 
ing observing mode was selected: field of view of 3x3"^ and 
spaxels of 0.100x0.050"^, wavelength range from 1.95 /im to 
2.45 fim and spectral resolution of 4000, corresponding to an 
instrumental broadening of FWHM=5.5A (this translates into 
FWHM=75 km s"' or cr=30kms"'). Since the angular distance 
between the clusters is about 4", we had to point on each tar- 
get separately. The total integration time for each pointing was 
set to 1800 s. The observations were carried out under regu- 
lar and stable atmospheric conditions, with a seeing of 1.2", a 
Strehl ratio of 42%, leading to an angular resolution of about 
0.2" {A/D-Q.Q5" at 2.1 fim), measured on the data, as detailed in 
Sec. [22] 

An attempt to fully reduce the data cubes using the ESO 
pipeline (Modigliani et al. 2007) was first performed, but left too 
strong sky line residuals on the final spectra. These residuals re- 
sult from the change of sky line intensity between the observa- 
tion of the object (OBS) and that of the sky (SKY). P-Cygni-like 
profile residuals are also observed in some cases. They are due to 
a slight shift in the wavelength calibration between the OBJ and 
S KY telescope p ositions. These effects are thoroughly discussed 
in iDaviesI (|2007|) and specific routines have been developed to 
correct them. They were kindly provided to us by R. Davies and 
allowed to improve the sky subtraction. All details about these 
procedures can be found in .Davies. (20071) . 

We used the following tools in the reduction process: the 
ESO pipeline, the routines by R. Davies, our own IDL rou- 
tines, and the data visualisation tool QFitsView, developed by 
Thomas Ott at MP^ 

The observation of each target was split into two blocks (two 
OBs), each of which is a sequence of six integrations alternat- 
ing object and sky (OBJ-SKY-SKY-OBJ-OBJ-SKY). For each 
field, we hence had six pairs (OBJ+SKY) of cubes. A standard 
jittering procedure was used for each sequence. 

The data reduction was performed as follows. First, the raw 
frames were cleaned from the bad lines that are created by the 
data processing hardcoded at the detector level. The correction 
was performed using the routine provided in the SINFONI data 
reduction cookbook. We then reduced the calibration frames us- 
ing the pipeline. Second, pipeline basic reduction was applied to 
each individual science cube (and not to each individual pair of 
cubes, as would be the case for standard pipeline reduction): the 
data were ffat-fielded and corrected for dead/hot pixels; interpo- 
lations to linear wavelength and spatial scales were applied, af- 
ter which the slitlets were aligned and stacked to create the final 
cube. The reduced sky cubes were then subtract ed from t he cor - 
responding object cube using the routines from IDaviesI ( l2007l) . 
A sample of a spectrum without/with residual correction is dis- 
played in Fig.[T] It shows that we did not entirely get rid of the 
artefact P-cygni like residuals. Hence, for enhancing the legibil- 
ity of the data, we decided to perform a cosmetic step after the 
correction, that replaced the pixel values by the median value 

' program ID: 076.B-0372(A) 

" http ://w w w. mpe. mpg . de/~ ott/QFi ts View 
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over the wavelength range of the main lines. These ranges are 
highlighted in grey on the spectra displayed in the paper. The 
only line of interest that is critically affected is H2 3-2 S(3). 




2.17 2.18 2.19 2.20 2.21 2.22 

wavelength [/^m] 
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Fig. 1. Example of the sky-subtraction correction, for an average 
case. We display a portion of the 1" radius aperture spectrum of 
one of the M4 data cube pairs. On top, the P-Cygni-like residuals 
are stronger before correction, and below, after correction, the 
residuals have been attenuated. The emission line is Bry. 



Then, for each pointing, the six sky-subtracted reconstructed 
cubes were combined. For this, relative astrometry was per- 
formed between the cubes, using the header information about 
the jitter sequences. 

We simultaneously performed the relative wavelength re- 
sponse and telluric feature corrections using the ESO provided 
data products of 27 observations of telluric standard stars. A me- 
dian response, corrected from the stellar feat ures, wa s then com- 
puted using the stellar spectra published bv iPicklell (Il998l) that 
are available on the ESO website. A scaling factor, following 
the Iraf telluric routine, was applied to adjust the depth of 
the atmospheric features. We also used this set of standard star 
data to evaluate the uncertainty on the flux calibration, which 
we call hereafter "photometric uncertainty", in contrast to "mea- 
surement uncertainty". The standard deviation of the conversion 
factors for these star flux measurements is 40%. Because these 
observations were obtained on different nights, we can be confi- 
dent that the "photometric calibration uncertainty", due to both 
the uncertainty on the standard star fluxes and the sky trans- 
parency variations, is at most +20%. 

To perform the relative registration of the cubes, we first 
built K-band images by convolving by the filter transmission, 
and then degraded the image and pixel size to adjust it to that of 
the ISAAC K-band image published by GA08. For each image, 
we ascertained the world coordinates of the brightest peak using 
the ISAAC image to register the cubes. We hence used the same 
astrometric reference as in GA08, where the [0";0"] coordinate 
was set to the peak of the AGN in the VISIR 12.8 ;um image. 

Once the registration was performed, we use the comparison 
between the degraded SINFONI K-band reconstructed images 
and the ISAAC images to perform the absolute flux calibration 
of the images. This choice implies that the flux measurements 
on the ISAAC and on the SINFONI data agree by definition. 

The last step of the reduction is the correction of the wave- 
length calibration to set the radial velocities in the local standard 
of rest (LSR). 



3. Data analysis 

This section presents the methods used to extract the quantitative 
information from the dataset. 

3. 1 . Map reconstruction 

The first step of the analysis of the reduced data cubes is re- 
constructing 2-D maps from 3-D cubes. For three chosen emis- 
sion lines (He I 2.06 /vm, H2 2.12/im and Bry 2.16;um), we per- 
formed an automatic Gaussian fit for each spectrum of the cubes, 
and built 1) the radial velocity maps (Gaussian profile central 
wavelength), 2) the radial velocity dispersion maps (cr of the 
Gaussian), and 3) the Gaussian flux maps from the fit parameter 
values. These maps are displayed in Fig. [S] from left to right, 
line flux, line radial velocity, and line width. A K-band map is 
obtained by convolving each spectrum by the filter transmission 
curve. A narrow band (excluding any spectral feature) map is 
also computed; it shows that the difference to the full K-band 
map is negligible. Accordingly, we only considered the full K- 
band maps. The K-band reconstructed map is displayed in Fig.|3] 
in the top left frame. This figure will be described in Sec. 13.21 
when we define the set of apertures used for performing the mea- 
surements. 

The Bry large-scale kinematics (radial velocities) are consis- 
tent with the CO kinematics. The Bry velocities can be compared 
to the large-scale isovelocity contours in SA07, which clearly 
show a rotation pattern with non-circular motions. From their 
tables, the measured LSR velocities (at the peak of brightness 
temperature for '^CO(2-l)) are 1560kms"', 1540kms"'and 
1480 km s"' for the CO knots associated with M4, M5 and 
M6. In the SINFONI data, for a 1.0" diameter aperture, we 
find 1555kms-i, 1520km s^^and 1475 km s"' for Bry, and 
1560 kms"', 1530 kms"', and 1470kms"' for H2 respectively. 
We measured these velocities with a simple Gaussian fit at the 
positions [0.41;7.12] , [2.80;9.86] , [4.80;6.90]. These results in- 
dicate that the Bry velocities are slightly blueshifted with re- 
spect to CO. The '2CO(2-l) FWHM are 70 km s"', 60 km s"' and 
80 km s~^ respectively, while Bry is broader with an FWHM of 
145 km s"', 145 km s"' and 200 km s"' respectively. 

3.2. Measurement procedures 

To perform measurements on the data cubes, we chose to first 
define a set of significant apertures to reduce the amount of in- 
formation to be analysed. We used the K-band, Bry, and H2 re- 
constructed maps to select the positions of these apertures. We 
chose the most prominent sources in the three images, as well 
as a few other apertures that sample regions of lower surface 
brightness. These apertures are drawn and labelled over the re- 
constructed maps in Fig.|3] which show from left to right, the K- 
band, Bry and H2 maps. Below each map we show the set of la- 
bels used for identifying the apertures: a number increasing with 
Right Ascension. Their radius of 0.1" allows one to sample sin- 
gle angular resolution elements. The measurements performed 
are consequently the least possibly affected by source mixing. 
Results for most measurements described below are given in ta- 
ble [1] |2] |3] and |4] and can be used for subsequent analysis. The 
line fluxes given in the tables have not been corrected for any 
extinction. 

The uncertainties on the different measurements have two 
origins: the uncertainty due to the noise and the measurement 
procedure, and the uncertainty on the calibration of the data, 
which we called photometric uncertainty in Sec. |2] In the ta- 
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H„ tlux [10"" erg s"' cm"'' arcsec"''] H„ Vrad [km s"'] H„ dVrad [km s"'] 
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Fig. 2. Gas emission around the extremely massive clusters in NGC1365, observed with VLT/SINFONI. First row: Bry flux, radial 
velocity, and velocity dispersion, from left to right. Second row, same for H2; third row, same for He I. The axes are in arcsec, north 
is up and east is left. 



bles, we only give the measurement uncertainties. The total max- 
imum uncertainties for the fluxes can be obtained by quadrati- 
cally adding the +20% uncertainty derived in Sec.|2]to the mea- 
surement uncertainty. 

To measure the angular resolution on each cube, we chose 
sources that are bright and narrow and measured their profiles 
on the K-band reconstructed image. We chose four sources: a 
telluric star, the K-band bright cluster Nl (to the west of M4), 
M5 and N2, the bright cluster on the southern part of the M6 
field. The profile measurements are shown in Fig. |4] The cuts 
are obtained along PA=0°, 45°, 90° and 135°, at the position 
of the source given in table [T] and at the centre of light for the 
star Since the data are AO-corrected, it is natural for the tel- 
luric star to have a narrower profile, because the Strehl ratio 



is higher In the best case, the achieved resolution would hence 
be 0.1" FWHM. For our dataset, the values measured are 0.23", 
0.20" 0.25", for the M4, M5 and M6 cubes, respectively, and 
are to be considered as upper limits, since we cannot be certain 
that the sources are not resolved. We used the fits of the spatial 
profiles for comparison with the profiles measured on the data. 
In Fig.|5] we compare the profiles of the main sources (M4, M5 
and M6) for K-band, Bry and H2 2-1 S(l) to the measured refer- 
ence profiles to check if they are resolved. 

The centre of each aperture was first chosen by eye on the re- 
constructed maps, and then taken as the centre of light in a box 
of 3x3 pixels around the initially chosen position. For apertures 
that coincide with emission peaks, the given aperture position is 
a measurement of the source position. For apertures defined on 
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6420 6420 6420 

KA [arcsec] HA [arcsec] HA [arcsec] 

Fig. 3. Top row, left to right: reconstructed maps for K-band flux density, Bry flux and H2 1-0 S(l) flux. The chosen apertures are 
drawn on the maps. The aperture labels are given for each map on the empty map below it for clarity. The label number is given in 
increasing Right Ascension. North is up and east is left, and the X and Y axes are in arcsec. 



extended emission regions, the precise position is not meaning- 
ful. Different peak positions coming from measurements inside a 
single cube (for example the K-band peak and a Bry peak) have 
a good relative registration precision, up to a fraction of pixel 
size (1/5 of a pixel corresponds to 0.01")- Between cubes, the 
registration was performed using the ISAAC images, and hence 
has the precision of a fraction of an ISAAC pixel size (1/5 cor- 
responds to 0.03")- 

The emission line fluxes through the given apertures are 
derived from automatic Gaussian+continuum fitting of the 
line profile s , carr ied out with the idl routine mpfitpeak by 
lMarkward3 (l2009l) . We measured the following lines: BiS, Bry, 
He I, and H2 roto-vibrational 1-0S(3), 1-0S(2), 2-1 S(3), 1- 
0S(1), 2-1 S(2), 1-0 S(0), and 2-1 S(l). The errors on the fluxes 
and widths were derived from the routine results, whereas the 
error on the line centre was simply taken to be 1/5 of the line 
FWHM. For the flux measurements for lines below 2 jum, the un- 
certainty is greater due to the uncertainty on the response, be- 
cause of the poorer atmospheric transmission in this region. We 
evaluated the uncertainty using the spectrum of M6 (apertures 
31 and 32), which displays a steep continuum. If we extrapo- 
late the continuum slope down to 1.95 yum, we find 0.0015 mJy, 
while we measure 0.0025 mJy on the spectrum. This means that 
the error can be higher than 50%. 

We measured the equivalent width of the CO 2.29 fim ab- 
sorption band (EW(CO)) using the following procedure. First, 
we built a spectrum of the background emission of the region. 
Using a histogram of the K-band flux density, we assumed that 
the pixels with less than 0.7mJy arcsec"^ belong to the back- 
ground. A median background spectrum was then built for each 



of the three cubes. For each aperture, we measured EW(CO) on 
the background-subtracted spectrum. The continuum fit was per- 
formed over three intervals between the bands, while the depth 
of the feature is measured between 23050 A and 23 160 A. The 
measurement error given in table [T]corresponds to + 0.3cr of the 
continuum around the fit. If the continuum signal to noise ratio 
(S/N) was less than 3, we only considered the upper limit of the 
EW(CO) (feature width divided by 0.3cr). 

At the defined aperture positions, we performed larger aper- 
ture background-subtracted flux measurements (aperture radius 
of 0.2" and background in an annulus between 0.2 and 0.6"). 
We used the djs_phot for the measurement and error evalua- 
tion. From these two measurements, we extract what we refer to 
in the following as the equivalent width of Bry: EW(Bry). 

We measured the slope of the continuum (after background 
subtraction) as the ratio between the fluxes measured through 
two artificial rectangle 200 A broad filters, which do not include 
any spectral feature and are centred at 20300 A and 22100 A. For 
these flux densities, the error was set to half the dispersion of the 
continuum. 



4. Observational results 

This section describes the main observational results derived 
from the analysis of the data described above. These results are 
discussed in Sec. |5] First, we present the results related to the 
images reconstructed from the data cubes and then those derived 
from the spectral measurements. 
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R flux [10 ^ mjy arcsec '] K flux [mjy arcsec '] B7 flux [10 erg s ' cm ^ arcsec ^] H,, flux [10 erg s ' cm ^ arcsec ^] 
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Fig. 6. Zooms on the 1.2x1.2"^ around M4 (first row), M5 (second row) and M6 (third row). For each source, from left to right we 
display the R-band image (HST), the K-band image emission, the Bry emission and the H2 l-OS(l) emission. The crosses mark the 
position of the Bry peak. North is up and east is left, and the X and Y axis are in arcsec. 



4.1. Images 

Fig.|6]shows the zoomed 1.2x1.2"^ fields around M4 (top row), 
M5 (middle row) and M6 (bottom row). From left to right we 
display (1) the HST "R-band" image (filter F606W), (2) the 
SINFONI K-band image, (3) the Bry image and (4) the H2 1- 
0S(1) image. The crosses mark the position of the Bry peak as 
defined by the procedure described in Sec. 13.21 It is important 
to keep in mind that the pixels on the cubes before combination 
are rectangular spaxels, which appear as 0.05x0.1"^ rectangu- 
lar vertical pixels on the reconstructed images. This can give the 
illusion that some sources are extended in the north-south direc- 
tion. 

The cluster positions were defined in GA08 as the peak of 
[Nell] 12.8 yum emission. Now, we can refine this positioning 
by taking the Bry peak position on the 0.2" resolution images. 
The [Ne II] and the Bry positions are consistent, and the differ- 
ence for M5 arises because, at the resolution of SINFONI the 
source is resolved into two bright knots. The new positions for 
the clusters are [0.39";7.13"] for M4/source 10, [2.80";9.80"] 
for M5/source 18 and [4.78";6.93"] for M6/source 31. The ab- 



solute precision of this registration with respect to the AGN is 
0.03" (+0.015", see Sec, fell. 

At the new achieved angular resolution, we see that the 
Bry and H2 sources that seemed coincident in the former ob- 
servations at lower resolution are in fact spatially distinct emis- 
sion knots. It is expected that the formerly superposed Bry and 
H2 knots are not anymore perfectly coincident when the angu- 
lar resolution is sufficiently high. This has previously been ob - 
served with a similar technique in IIZw40 bv lVanzi et al.l (l2008h . 
In our case, the data show that M4 is marginally resolved in K 
and Bry, while the associated H2 peak is resolved (FWH]yi~0.3") 
and is shifted by about 0.1" to W-SW. A secondary H2 source 
(source 3) is detected about 0.5" to the east of M4, and has no 
counterpart in any other emission. M5 is clearly resolved into 
four sources (apertures 16, 18, 21, 22) following a cloverleaf 
configuration inside a circle of about 0.5"diameter Apertures 
18 and 21 correspond to an emission peak in K and Bry. These 
sources are not resolved individually. The associated H2 emis- 
sion lies outside this cloverleaf, in two resolved knots aligned 
along PA~70° and at a distance of about 0.8" (apertures 23 and 
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Fig. 4. Light profiles (ID cuts through the images) for a telluric 
star, Nl, M5 and N2. For the first three, we took cuts along 
PA=0°,45°,90° and 135°(black crosses, cyan asterisk, purple tri- 
angle and brown diamond). The median profile is drawn with a 
black line, and the Lorentzian fit with a salmon line. For N2, we 
only have PA=0° since it is on the edge of the image. The fit 
FWHM are 0.10", 0.23", 0.20" 0.25". 
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Fig. 5. Profile cuts for M4, M5 and M6 (position given is the 
table for the light centre on the K-band image) along PA=45°, 
60°, and 1 10°. The dotted line is the comparison profile for each 
cube (as explained in the corresponding paragraph). It has been 
rescaled to give a meaningful comparison. 



15). M6 is detected and resolved in K and Bry. It has an asym- 
metrical shape, with an extension towards the S-E, which can be 
seen as an extended wing in the profiles in Fig.|4] The source as 
a whole is detected in H2, but the maximum of H2 is displaced 
by about 0.1" to the south with respect to the Bry emission. 

Fig. 12] shows the full combined gas emission reconstructed 
maps for Bry (top row), H2 1-0 S(l) (middle row) and He I (bot- 
tom row). From left to right we display the line flux, line centre 
and cr maps. We detect extended emission on scales of several 
arcsec, both in Bry and H2. They have the same global distribu- 
tion, but differ at the sub-arcsec scale: M4 and M5 are located 
on a filament-like emitting structure oriented along PA~45° and 
extending over about 7", with a width of about 1". A distinct and 
more symmetrical emission region of about 3" diameter is ob- 
served around M6. The bright emission peaks are localised well 
in the middle of the extended emission regions. Apart from M4, 
M5 and M6, and their closely associated sources, there are only 
few secondary weak peaks (apertures 4, 7, 13 and 28 for Bry). 

4.2. Spectra 

Fig. I2] shows several full-wavelength range spectra of (a) the 
background as defined in Sect. 13.21 (b and c) the two brightest 
K-band clusters, and (d, e and f) M4, M5 and M6. Over the three 
first spectra, we plotted in grey norma l star spectra taken from 
flie library of IWaUace & Hinkld (Il997h for comparison. The line 
wavelengths given in the following are the rest wavelengths, the 
observed redshifted ones are about 0.01 fim longer. 

The full-field combined K-band map is displayed on the left 
of Fig. |3] What we defined as background continuum in Sect. 13.21 



appears in red (up to 0.7 mJy arcsec"^). The corresponding spec- 
trum is given in Fig. |7]a. Across this background, we mostly 
detect continuum sources along the border of the dust lane, on 
its western side. Only a few sources are detected on the dust lane 
itself, except for M6: we identified the sources corresponding to 
apertures 5, 12, 24 and 25. Two bright K-band clusters are ob- 
served: apertures and 26. Their spectra are given in Fig.|7j5,c. 

The emission lines initially present on the background spec- 
trum were removed artificially, and a pure continuum spectrum 
is left. In grey we overplot the spectrum of HR6705 (K5III). 
The next two spectra (apertures and 26) correspond to the two 
brightest K-band sources. From these spectra, no emission lines 
were removed. The overplotted stellar templates are HR2508 
(Ml -1- Ib-IIa). The spectral slopes of the templates were slightly 
reddened to adjust them to the observed ones. The spectra dis- 
play strong CO absorption as well as metallic line absorptions. 
The next three spectra are M4, M5 and M6 (apertures 9, 18 and 
31). The three sources show a rising continuum towards the red 
and several emission lines, labelled on the spectrum of M6: Bt6, 
Bry, the 2.06 He I yum line and seven H2 roto-vibrational lines. 
For M5 and M6, a broad He I 2.11 yum observed close to the 
blue side of H2 1-0 S(l) is detected. The Gaussian fit FWHM is 
50 A and 38 A for M5 and M6, respectively, while it is about 
10 A for the other lines. Weak He II absorption at 2. 19yum is de- 
tected for M5 and M6. For M6 only, we detect two weak [Fe III] 
emission lines at 2.22 jum and 2.24 yum. An unidentified emission 
line is detected on the spectrum of M6 at an observed wavelength 
of 22097 A, which corresponds to a rest wavelength of 2 1 990 A. 
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Fig. 7. Top to bottom: (a) spectra of the background (defined in Sec. 13.21 1 + HR6705 in grey, (b) aperture + HR2508 in grey, 
(c) aperture 26 + HR2508 in grey, and (d) aperture 9 (M4), (e) aperture 18 (M5) and (f) aperture 31 (M6). (d), (e) and (f) are not 
background subtracted. 



For M4, M5 and M6, we detect an asymmetry in the Bry pro- 
file. The same asymmetry is observed for He I 2.06 ;um. Fig. [8] 
displays the profiles for M4, M5 and M6 in the left column from 
top to bottom, and for a set of off-centred apertures (4, 13 and 
28) in the right column. The intensities are shown on a logarith- 
mic scale. The null velocity is computed as the central velocity 
of the Gaussian fit of the line. For the positive velocities, the 
crosses just show the data, and these are mirrored on the side of 
the negative velocities, from which one can appreciate the de- 
gree of asymmetry of the profile. Clearly, M4, M5 and M6 have 



an extended blue wing, while the spectra from the three other 
apertures display symmetrical profiles. 

Fig. |9] gathers a set of diagrams that plot various observables 
measured on the data following the procedures detailed in the 
former section. For all diagrams, the blue dots mark the bright 
Bry sources, and the red dots mark the H2 emission knots. The 
purple labels refer to either the aperture numbers as defined in 
Fig.|3]or directly the source names for M4, M5 and M6. As de- 
tailed in Sec. 13.21 the equivalent width of the CO 2.29 /^m band 
was measured through all defined apertures. The EW(CO) vs. 
the aperture number is plotted in Fig.|9]a. Absorption is detected 
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through all apertures, except for the Bry peak sources M4, M5 
and M6, for which the upper limit is a few A. For the sources 
that are not bright in Bry, the EW(CO) is about 20±5 A. Fig.|9]b 
shows as a histogram the number of apertures per EW(CO) bin 
and identifies two groups. In Fig. |9]c, we plot the EW(CO) vs. 
the EW(Br7). We see that the two groups are well identified 
in both quantities, with a well-defined anti-correlation between 
EW(Br7) and EW(CO). M4-2 (or aperture 6, the K source right 
next to M4) has a high EW(CO), and its also high EW(Br7) is 
likely due to its being superimposed over the extended Bry emis- 
sion from M4. The EW(Bry) of the older sources range from 
OA to 20 A, and their EW(CO) from 15 A to 25 A. The other 
sources have higher EW(CO), ranging from 1 10 A to 150 A, and 
low EW(Bry). 

From the line analysis of our dataset, the only possibility to 
evaluate the extinction is to use the ratio between the two hy- 
drogen lines that we can measure: Br<5 and Bry. Unfortunately, 
as stated in Sec. 13.21 the uncertainty on the flux calibration in 
the region of Bid is greater than in the other regions of the 
spectral range. Indeed, if we blindly measure the Bid fluxes, 
we find the Bi6/Biy ratio to be higher than that of the theo- 
retical case B (0.66). We decided to apply a correction factor 
of 0.0015/0.0025=0.6 to the Bi6 fluxes (according to the contin- 
uum shape, as explained in Sec. 13. 2b . In Fig.|9]d, we plot the con- 
tinuum slope as defined in the previous section vs. the BiSfBiy. 
The correlation suggests that the origin of the continuum slope is 
the extinction. We transformed both quantities in Av in Fig.|9]e. 



To derive the extinction from the continuum slope, we consid- 
ered that the flatter source (N2 or source 26) has no extinction. 
We see that at least for the brightest Bry sources (cyan dots), 
the values are well consistent and are in an Av range from 10 to 
30 mag. 

Fig.|9]also displays three additional plots. Plot|9]f shows that 
there is a possible correlation between the width and flux of Bry. 
We will not discuss any implications of this because the apparent 
correlation is weak, but we think it is worth mentioning, given 
the number of results in the literature about the L-cr relation in 
dwarf galaxies (iBordalo & Tellesll20T 1, and references therein). 
Plot|9]g shows that the H2 and Bry line widths are correlated, 
with H2 being narrower than Bry (H2~ 0.85Bry). Interestingly, 
in the "filament" joining M4 to M5, the Bry FWHM values are 
about 10 A, while they exceed 14 A in the region of M6. This 
region of higher line width for both Bry and H2 also displays a 
blueshift higher than the surrounding, as can be seen in Fig. |2] 
The last plot (Fig. |9]h) shows that the H2 fine ratio 2-1 S(l)/1- 
S(l) does not depend on the aperture through which it is mea- 
sured and has a value of ~0.2+0.1. This is the ratio of the fluxes 
without any correction for extinction. The effect of extinction 
on this ratio is small since the lines are in the NIR and their 
wavelengths are close. For an Av of 30 mag, which represents 
the most extreme possibility according to Fig. |9]e, an observed 
ratio of 0.2 implies an actual ratio of 0.15, within the error. 



5. Discussion 

In this section, we discuss the results derived from the dataset 
and presented above. The improvement of the spatial resolution 
with the SINFONI data allows us to re-examine the conclusions 
from GA08, some of which have been summarised in the in- 
troduction. First, we show that the new data allow us, contrary 
to GA05 and GA08, to resolve the K-band stellar continuum 
emission, and hence prove that the sources are no single clus- 
ters but rather cluster complexes. Then, we discuss the implica- 
tions of the non-detection of the CO absorption bands, contrary 
to our previous work, derive new ages for the clusters, and dis- 
cuss the apparent inconsistency with our previous estimate of 
7 Myr, based on the MIR line ratios. In Sec. 15.21 we discuss the 
origin of the extended gas emission around the clusters, making 
full use of the integral field capacity of SINFONI. 

5.1. Multiplicity and age 




age [Myr] age [Myr] 

Fig. 10. Starburst99 equivalent widths for Bry emission on the 
left, and CO 2.29 yum absorption band on the right. 
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First, the SINFONI data show that M4, M5 and M6, which 
we repeatedly referred to as clusters, are resolved in the K-band 
map at a resolution of 0.2", which corresponds to about 20 pc 
(Sec. 14. H and Fig. ID. Since the K-band emission traces the stel- 
lar photospheres, these sources are very likely cluster complexes. 
For M4 and M6, we do not detect multiple sources, but resolve 
the peaks, hence we cannot formally refute the possibility that 
they are single large "clusters". We note that the K-band source 
M4-2, which at first sight could be considered to belong to the 
M4 cluster complex, is indeed an older source, with deep CO 
absorption (see next paragraph). For M5, we detected four un- 
resolved sources in an area of about 0.5x0.5"^=45x45pc^. The 
observational evidence that M4, M5 and M6 are cluster com- 
plexes is the crucial new piece of information that the SINFONI 
data reveal. 

A second important piece of information is gained because 
the SINFONI data allow us to reliably probe the presence of 
the CO bands at high angular resolution, after careful evaluation 
of the background emission. Since the rapid evolution of very 
massive stars turns them into red supergiants in a few million 
years, the NIR continuum is dominated by them from 6-7 Myr 
onwards and the deep CO bands appear. Fig. [TO]shows the evo- 
lution predictions of EW(CO) between and 10 Myr for a single, 
normal IMF, Z=0.04, instantaneously formed cluster. It shows a 
st eep rise between 6 and 7 Myr. As shown by Fig. 101 and 102 
of iLeitherer et al. l (fT999l) . the beginning of the steep rise of the 
CO depth will not depend on the precise values of these parame- 



ters, since it only depends on the evolution of very massive stars. 
For M4, M5 and M6, the measured upper limits on EW(CO) 
of 2 to 4 A (Sec. 13.21 and Fig. |9]a) formally imply an age up- 
per limit of 6.5 Myr. For M5-2 and M5-4, we also have a lower 
limit for EW(CO) and the corresponding ages lie between 6 Myr 
and 7 Myr. According to Starburst99, EW(CO) should never ex- 
ceed 17 A. Fig. |9^,b show higher values, but given the order of 
magnitude of the measurement error, this is not significant. The 
EW(CO) of the two brightest K-band sources (apertures and 
26) lie somewhere between 19 and 23 A, and hence must be 
at the maximum of EW(CO), corresponding to an age between 
8 and 30 Myr for instantaneous formation, and between about 
10 Myr and 100 Myr for continuous formation. 

Another common age indicator for very young clusters is 
the equivalent width of Bry, EW{BTy). This is valid under the 
hypothesis that we are able to actually measure the continuum 
from the ionising cluster only, and the Bry flux emitted by the 
gas excited by this cluster only. Moreover, we must assume that 
the extinction towards the stars and the emitting gas is the same, 
which is not obvious since the spatial scale of the emitting gas 
is expected to be larger than that of the stars. For this last re- 
mark, we can object that this does not seem to be the case for 
the small apertures for the Bry-emitting sources (M4, the com- 
ponents of M5 and M6) since we measure a very good corre- 
lation between the extinction derived from the K-band contin- 
uum slope and that derived from the Bi6/Biy line ratio (Sec. 14.21 
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and Fig.|9]e): the cyan dots lie close to the dotted line. Fig.|9]c 
plots the EW(CO) vs. the EW(Bi-y) (measured following [3!2] l 
and shows two groups, the older sources, with low EW(Btj) and 
high EW(CO), and the younger sources with high EW(Bry) and 
low EW(CO). For M4 and M6, EW(Brr) is about 130 A, and 
it is about 115 A for the M5 sources. The measurement errors 
are small because of the measurement procedure, but, given the 
uncertainties mentioned above, the real uncertainty for the quan- 
tity that it makes sense to compare to model predictions must be 
larger Nevertheless, even if we consider a large error of 50 A, the 
formal ages implied by comparison with the Starburst99 models 
(Fig. [Tol l lie in the narrow range from 4.5 Myr to 5 Myr 

In this paragraph, we briefly state why other potentially 
powerful age indicators cannot be used in our particular case. 
(1) Wolf-Rayet features: We can easily see that WR features 
in the spectra of the clusters, if present, would be too faint to 
be detected. To estimate the flux of a WR in the K-b and, we 
used the WR 102ca spectrum of iHomeier et al. I (l200l . which 
is is at a distance of 8.5 kpc (quintuplet cluster) and has a flux 
of 0.04Jy. According to Stai-burst99, the NGC1365 7Myr old 
clusters should contain about 1000 WR stars. At the distance 
of NGC1365 (18Mpc), this population would have a flux of 
only O.OlmJy, undetectable in our data. (2) Broad-band colours: 
NGC1365 has been extensively observed by HST, and hence im- 
ages in different filters are available from the archive, and could 
be used to measure the broad-band colours of the clusters and 
derive ages. M4 and M6 are heavily obscured, and hence not de- 
tected in the visible, but M5 is detected. Unfortunately, GALEV 
models (Kotulla et al. 2009) show that we would need a O.lmag 
precision on the Av of the cluster itself to distinguish between 
4Myr and SMyr clusters. (3) He I to Bry ratio: In principle, this 
ratio would give information about the hardness of t he UV con- 
tinuum emitted by the cluster, and link it to its age (Dovon et al.' 
[ 1992), but this has been demonstr ated to be unreliable by Shields 
(1 19931) and lLumsden et alJ (l2003h . 

In summary, for a standard instantaneously formed 
Starburst99 model, there is inconsistency between the different 
mass estimates: 

1 . The non-detection of the [Arlll] and [SIV] lines shows that 
the most massive, hence hard UV emitting stars, have died 
already, hence that the cluster is older than 7 Myr (Fig. 15 
and Fig. 16ofGA08). 

2. The non-detection of the CO absorption bands shows the ab- 
sence of red super giant stars, hence proves that the clusters 
are younger than 6 Myr. 

3. The EW(Bry), even though we believe is not the most reli- 
able age indicator in our case, is consistent with ages in the 
range 4.5 Myr to 5 Myr 

Choosing a more complex star formation history cannot 
solve the above contradiction since the age limits measured 
by the non-detection of the MIR lines and CO lines are in- 
deed the time that has elapsed after the last massive star for- 
mation episode ended. We trust our previous MIR observations, 
since we have an example of a similar or even slightly lower 
quality spectrum for a comparable, but younger object, an d the 
rSIV] line is much brighter than [Nell] jMartm-Hernandez et aTl 
l2005h . as detailed in the introduction, and predicted in GA08. 
The VISIR spectra were of good quality, several expected spec- 
tral features are detected, and the [Nell] line has a high S/N ratio. 

An easy way of making the Starburst99 model consistent 
with the observational constraints is to lower the upper mass cut- 
off" of the IMF (Mmax) to values between 25 Mq to 30 Mq. In this 



case, the cluster UV continuum is simply never sufficiently hard 
to maintain the ionisation of [Arlll] and [SIV]. 

For Mtnax=25Mo and Minax=30Myr, we find that the ob- 
servational constraints are consistent with age ranges of [6.2- 
6.4] Myr and [5.62-5.64] Myr, respectively. Hence, under the hy- 
pothesis of a "low" upper mass cutoff for the IMF, the ages 
derived might be only slightly lower than our first estimate in 
GA08, in the range [5.5-6.5] Myr. The corresponding mass esti- 
mates for the clusters would then be between 3 and 10x10^ M©. 
We do not claim any evidence for a low upper mass cutoff IMF, 
but only show that the observational constraints are consistent 
with it, and inconsistent with the "standard" Starburst99 model. 
Given the large masses of the these systems, the most massive 
stars cannot be absent just on the basis of stochastics and an 
astrophysical mechanism should be invoked to explain the pos- 
sibility of an actual absence of stars with masses greater than 
25 Mq. This is beyond the scope of the current analysis. 

In addition to the absolute ages of the cluster complexes, we 
note the presence of spectral features that might witness age dif- 
ferences between the clusters. We see that the slope (and extinc- 
tion) increases between M5, M4 and M6. Interestingly, in the 
spectrum of M6, we detected weak [Fe III] lines that may reveal 
a slightly harder ionising field, hence a slightly younger age. The 
younger cluster complex being the most extincted is consistent 
with the idea of ongoing gas expulsion at that stage of evolution. 

Even though we cannot derive any precise extinction from 
the SINFONI cubes, we estimated of about 10 to 30 mag through 
the small apertures extinctions, with a good consistency between 
the extinction from the g as emission lines and the star continuum 
(Sec. EJand Fig|9](e) ).'Ganiano et al.'(2008') derived compara- 
tively low extinctions towards M4, M5 and M6 of 13.5, 3.2 and 
8.5 mag, through much broader apertures. 

In summary, the new dataset has highlighted that M4, M5 
and M6 are cluster complexes, and provided evidence (EW(CO) 
and EW(Bry) at high resolution) for ages (5.5 Myr to 6.5 Myr) 
that are slightly younger than previously derived (7 Myr). We 
showed that the inconsistency between the ages derived from 
the EW(CO) and EW(Br7') on one side, and the MIR line ratios 
on the other, can be overcome by assuming an upper mass cutoff 
of 20 to 25 Mq for the IMF. We also showed that there is some 
indication for an age gradient across the sources, M6 being the 
youngest. The comparison between the old and the new extinc- 
tion measurements strongly suggests that the extinction is patchy 
and higher towards the main sources. 

5.2. The extended gas emission 

In this subsection, we discuss the origin of the extended gas 
emission around the cluster complexes. A direct qualitative com- 
parison between the maps for M4, M5 and M6 and the maps 
availa ble for the we ll-known giant H ii region 30 Doradus in the 
LMC dlndebetouw et al. 2009; Peflegrini et al. 2010, 2011) can 
help in understanding the origin of the extended emission in the 
observed field. An impressive image of the whole Tarantula neb- 
ula was published in the ESO Photo-Release eso0216. It shows a 
nebula more than 300 pc across, which corresponds to an angu- 
lar size of more than 3" at the distance of NGC1365. The ioni- 
sation of the 30 Doradus nebula is dominated by photoionisation 
from the central cluster R136. This shows that even though the 
presence of more uniformly distributed lower mass and recently 
formed clusters in the region of extended gas emission cannot 
be excluded, such a configuration is not necessary to explain the 
extent of the emission. The region in which we observe M4, M5 
and M6 possibly also only contains three giant HII regions, each 
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dominated by a central cluster complex confined in a ~ 50 pc 
region. That the emission peaks are located well at the centre 
of the extended emission in the case of M6, and close to the 
symmetry axis of the "filament" joining M4 to M5, strengthens 
the idea that the stellar content is dominated by central cluster 
complexes. The observed slight shifts in position between the 
emission peaks (continuum, Bry and H2) are indeed expected 
for a "30 Doradus"-like complex and clumpy gas structure for 
the following reason: at a resolution of about the total angu- 
lar size of the emitting region, these emissions are coincident, 
but the stars (sources of continuum), the ionised gas (source of 
Bry) and the warm gas (source of H2) occupy different regions 
of space, and may hence not be spatially coincident. This can be 
seen in the three colour image of R136 displayed in GA08 (Fig. 
12), with the V-band emission tracing the stars. Ha the ionised 
gas, and 8 fim the warm dust. When observed at a sufficiently 
high resolution, the difference in position between the peaks of 
these components can be unveiled. 

Indeed, the extended H2 emission is well consistent with 
the H II emission on scales of several arcsec, while they do not 
coincide at a higher level of details. Fig. [3] clearly shows that 
the H2l-0S(1) / Bry ratio increases with distance to the cen- 
tral sources. We may wonder whether this is a sign of local 
excitation of the H2 emission. The K-band spectra show a se- 
ries of H2 lines both for para- and ortho-hydrogen. The excita- 
tion mechanism for these lines can give clues about the pres- 
ence of an underlying population of smaller mass clusters. The 
low 2-1 S(l)/l-0 S(l)~0.2 ratio observed (Fig.|9]h) usually traces 
thermally excited H2 in shocks. However, for high densities of 
about 10^ cm"^ and for an intense incident FUV radiation field, 
coUisional de-excitation of UV- pumped rnolecules thermalises 
the lower energy lev els (Burton et al. ■ill 9901: f 



HoUenbach & Nattal 



119951: iDavies et all 12003,: ,Martm-HernandezetalJl2008h . That 
the central UV sources possibly maintain ionisation of hydro- 
gen at large distances suggests that the gas distribution is very 
clumpy, hence that the density is high. The only quantitative in- 
dication of the gas density comes from the observed ratio of the 
two forbidden [Felll] lines (^Gs H(, at rest 2.2183 /im and 
H4 at rest 2.2427 fim) detected in the spectra of M5 and 
M6 (apertures 18 and 31). We measured 0.38+0.05 for M6 and 
<0.5 for M5. If we compare thi s to the mode l predictions from 
iKeenan et all (Il992h reported in [Gilbert et"al] (12000). we find a 
density of 10^ cm"-'. This means that we cannot decide between 
the two excitation m echanisms from the 2-1 S(1)/1-0S(1) ratio 
only. As explained in lOavies et al. 1 12001 . while the v-2 vibra- 
tional levels are thermalised at high density, the v-3 level is ex- 
pected to retain a fluorescent population until even higher den- 
sities. In the transitions involving the v=3 level, the data only 
leave upper Hmits for 3-2 S(2), 3-2 S(4) and 3-2 S(5) while the 
3-2 S(3) line is blended with a skyline. For M6 (aperture 3 1), we 
found 1-0 S(l)/3-2 S(2) > 20, and for the "pure" H2 source (aper- 
ture 15), close to M5, we found a lower limit of 35. Considering 
a ratio 3-2 S(2)/3-2 S(l) equal to the standard para-ortho ratio 
of 1/3 (see next paragraph ), gives l-0S(l)/3-2S(l)>7 and 12, 
respectively. According to Martm-Hernandez et al. (l2008h . for 
densities of about lO^cm"^, the fluorescent ratio retains values 
lower than 8. The measured lower limits hence suggest that the 
H2 emission is thermally excited in shocks. 

As discussed by iMartm-Hernandez et al.l (l2008h . an ad- 
ditional diagnostic can be obtained by the measurement of 
the ortho-to-para ratio of the H2 molecules. Where shock 
excitation is dominant, this ratio <[), based on the stat isti- 
cal weights of the nuclear spins, is 3 dSmith et"an 1 19971) . It 
can take other values in the case of fluorescent excitation 



(ISternberg & Neufeldlll999h. F ollowing ISmith et al] (Il997h and 
iMartm-Hernandez et al.l (l2008l) . we can measure the ortho-to- 
para ratio for the v-1 states (Oi) measuring the fluxes of 1- 
0S(0), 1-0 S(l) and 1-0 S(2). The measured values of (Dj tend 
to be similar in all apertures, and slightly lower but consistent 
with 3, with a median of 2.4 an d an error of +0.7. As reported by 
Ma rtin-Hernandez et al.l (l2008l) . ratios in the range 1.5-2.2 have 
been measured in PDRs. Hence, the measured values are consis- 
tent with shock excitation, but have a median that could indicate 
part of fluorescent excitation. 

The different diagnostics for the excitation of H2 suggest a 
picture in which shocks are the dominant mechanism, while flu- 
orescence, probably also present, is marginal. That shock excita- 
tion at large distance from the central complexes is taking place 
might signal the presence of a lower mass population of young 
clusters. However, we can also predict the origin of these shocks 
in gas outflows produced by the central complexes themselves, 
without the need of invoking another "hidden" cluster popula- 
tion. The blue wing of the Bry line detected at the location of 
M4, M5 and M6 and displayed in Fig.[8]is possibly the evidence 
of such an outflow. The comparison with line profiles extracted 
away from the main sources, which are free of this asymmetry 
suggests that the outflows are indeed centred on the main sources 
M4, M5 and M6. The associated velocities are some hundreds of 
km s which means that given the size of the giant H 11 regions 
and their age, there has been enough time for gas heating through 
shocks to occur. The non-detection of the red wing can be simply 
attributed to extinction. 

In summary, we find that the new data are consistent with 
a picture in which the central cluster complexes constitute the 
source of excitation of the whole nebula through ionisation and 
shocks produced by gas outflows. Even though we do not discard 
it, we find no need to invoke a lower mass cluster population to 
interpret the data at hand. 



6. Conclusions 

Thanks to the integral field capacity of SINFONI, along with the 
improved angular resolution achieved with the adaptive optics, 
our knowledge and understanding of the three bright MlR/radio 
'embedded clusters' M4, M5 and M6 located in NGC1365 nu- 
clear region has greatly improved. 

First, the gain in angular resolution in the new K-band maps 
(compared to the former ISAAC K-band maps) allowed us to 
resolve the sources. This proves that they are not single clusters, 
but instead cluster complexes of several tens of parsec size. 

The possibility to resolve the spectral features at a high an- 
gular resolution allowed us to perform precise measurements of 
both the equivalent width of the 2.3 fim CO absorption bands 
and that of the Bry emission line, uncontaminated by back- 
ground or close-by surrounding emission. The unexpected non- 
detection of the CO bands together with the measured values of 
EW(Bry) suggest ages that are younger than previously derived. 
These new measurements, within the framework of a 'standard' 
Starburst99 cluster model, are incompatible with the MIR line 
ratios discussed in our previous paper. We do not have any re- 
liable explanation for this apparent inconsistency yet. But we 
note that it is a hint for an upper mass cutoff of the star clus- 
ter IMF lower than usually assumed with values below 30 Mq. 
In this case, the complexes would only be slightly younger (5.5 
to 6.5 Myr) than derived in our former work, with corresponding 
masses of [3-10]xlO* M© per complex. If we had kept a standard 
IMF, and simply discarded our previous estimate, we would have 



12 



Galliano et al.: cluster complexes in NGC1365 



found ages of 4 Myr, resulting in much lower masses, of about 
10^ Mq. 

The complete dataset available for these sources and the ex- 
tended emission around them suggests that this gas-rich region 
(SA07) consists of three co-eval multiple "30 Doradus"-type gi- 
ant H n regions, each powered by a central star cluster complex. 
There is a weak observational indication for a slight age gradi- 
ent among them, with the youngest (M6) displaying a steeper 
K-band slope and the presence of weak higher ionisation lines. 
A few secondary Bry maxima are detected, and witness the pres- 
ence of at least some lower mass clusters or cluster complexes, 
but we found no compelling need to invoke a fully populated 
power law cluster mass distribution to interpret the data at hand. 
The angular resolution achieved is sufficiently high to show that 
the distributions of the Bry and H2 line emission is not coinci- 
dent, indicating a patchy gas distribution, where the ionised gas 
is probably seen on the surface of denser knots, and reveals com- 
plex kinematics (line asymmetry and broad width). The colder 
H2 emitting gas is distributed all over the region, and appears 
to be predominantly thermally excited in shocks, rather than 
fluorescence in photo-dissociation regions. We interpret the ob- 
served extended blue wing in the Bry profiles as a signature of 
gas outflows from the main cluster complexes, and suggest that 
these outflows would be the natural source of shock excitation 
of molecular hydrogen over the entire region. 
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Table 1. Results from the measurements described in Sec. l3.2l for the apertures defined in the K-band image. The aperture number 
is defined in Fig. [3] and the corresponding name is given in the second column when appropriate. The positions of the apertures 
are (A(R.A.);A(dec)) in arcsec with respect to the position of the AGN in the MIR (see text). The K-band flux density is mJy. 
The Bry flux is in erg s"' cm"^. These two fluxes are background-subtracted fluxes measured on the reconstructed images, through 
"large" apertures as defined at the end of Sec. 13.21 (aperture radius 0.2", and background between 0.2" and 0.6"). The slope is 
defined in the last paragraph of Sec. 13.21 The equivalent width of the CO absorption band EW(CO) is measured through the small 
apertures following the procedure defined in Sec. l3.2l and is given in A. 
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Table 2. Measurements described in Sec. 13.21 for the apertures defined on the Bry image and the H2 image. The line fluxes are in 
erg s ' cm"^, and the line Gaussian cr are in A. 
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aperture 
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17 


1.7e- 


■15± 3.2e-17 


4.5± 0.2 


4.1± 0.1 


3.1+ 0.2 


4.1± 0.1 


4 




6.6e- 


16± 6.1e- 


■17 


4.9e- 


■16± 3.4e- 


17 


1.8e-16± 3.3e- 


17 


1.2e- 


■15± 3.2e-17 


3.7± 0.2 


4.1± 0.1 


5.3± 0.5 


3.9± 0.1 


7 




6.2e- 


■16± 7.1e- 


■17 


4.8e- 


■16± 3.6e- 


17 


1.5e-16± 1.9e- 


17 


l.Oe- 


■15± 3.1e-17 


3.7± 0.2 


4.7± 0.2 


2.3± 0.1 


4.0± 0.1 


8 




1.5e- 


■15± 8.9e- 


■17 


7.5e- 


■16± 4.4e- 


17 


3.8e-16± 3.7e- 


17 


2.4e- 


■15± 3.6e-17 


3.7± 0.1 


3.7± 0.1 


4.7± 0.2 


4.1± 0.1 


9 


M4 


1.3e- 


■15± 9.6e- 


■17 


7.3e- 


■16± 4.5e- 


17 


4.1e-16± 3.8e- 


17 


2.2e- 


■15± 3.5e-17 


3.8± 0.1 


3.7± 0.1 


5.0± 0.2 


4.2± 0.1 


11 




7.9e- 


■16± 8.6e- 


■17 


5.6e- 


■16± 3.8e- 


17 


3.0e-16± 4.8e- 


17 


1.3e- 


■15± 3.1e-17 


3.6± 0.2 


4.6± 0.2 


8.4+ 0.7 


4.0± 0.1 


13 





7.2e- 


■16± l.le- 


16 


5.1e- 


■16± l.Oe- 


16 


< 1.2e-15 




7.2e- 


■16± 5.4e-17 


5.4± 0.4 


9.2± 0.9 




4.9± 0.2 


15 




1.7e- 


■15± 9.4e- 


■17 


7.7e- 


■16± 5.2e- 


17 


2.7e-16± 3.2e- 


17 


2.3e- 


■15± 3.2e-17 


4.0± 0.1 


4.7± 0.2 


4.1± 0.2 


4.1± 0.1 


16 


M5-4 


l.Oe- 


■15± 1.2e- 


16 


4.5e- 


■16± 5.2e- 


17 


7.6e-17± 2.7e- 


17 


1.2e- 


■15± 4.8e-17 


3.5± 0.2 


4.7± 0.3 


2.3± 0.4 


3.7± 0.1 


18 


M5 


1.2e- 


■15± 2.1e- 


16 


7.4e- 


■16± 7.2e- 


17 


3.3e-16± 4.4e- 


17 


1.5e- 


■15± 5.5e-17 


4.1± 0.4 


5.2± 0.3 


3.5± 0.2 


3.4± 0.1 


21 


M5-2 


9.8e- 


16± l.le- 


16 


5.9e- 


16± 5.2e- 


17 


2.8e-16± 4.0e- 


17 


1.5e- 


■15± 4.5e-17 


3.1± 0.2 


4.5± 0.2 


3.6± 0.3 


3.7± 0.1 


22 


M5-3 


1.7e- 


■15± l.Oe- 


16 


7.4e- 


16± 4.9e- 


17 


2.4e-16± 3.6e- 


17 


2.1e- 


■15± 3.6e-17 


4.1± 0.1 


4.2± 0.1 


3.6± 0.3 


4.1± 0.1 


23 




1.8e- 


■15± 9.8e- 


■17 


7.1e- 


■16± 4.6e- 


17 


3.0e-16± 3.5e- 


17 


2.2e- 


■15± 3.2e-17 


4.1± 0.1 


4.1± 0.1 


3.5± 0.2 


4.2± 0.1 


28 




5.2e- 


16± 6.6e- 


■17 


3.2e- 


■16± 3.7e- 


17 


1.2e-16± 2.3e- 


17 


8.5e- 


■16± 3.2e-17 


4.8± 0.3 


3.7± 0.2 


2.4± 0.2 


4.1+0.1 


29 




l.Oe- 


■15± l.le- 


16 


7.4e- 


■16± 6.1e- 


17 


1.9e-16± 3.4e- 


17 


1.8e- 


■15± 5.1e-17 


5.1± 0.3 


6.9± 0.3 


4.7± 0.4 


5.0+ 0.1 


31 


M6 


8.0e- 


■16± 1.2e- 


16 


8.1e- 


■16± 8.0e- 


17 


1.8e-16± 4.8e- 


17 


1.6e- 


■15± 5.8e-17 


4.6± 0.4 


8.4± 0.4 


6.0± 0.8 


4.7+ 0.1 


33 




1.2e- 


■15± 1.4e- 


16 


7.1e- 


■16± 5.5e- 


17 


1.5e-16± 3.4e- 


17 


1.7e- 


■15± 4.3e-17 


6.4± 0.4 


7.0± 0.3 


3.4± 0.4 


5.7+ 0.1 


34 




5.0e- 


■16± 4.5e- 


■17 


3.5e- 


■16± 3.3e- 


17 


1.2e-16± 2.7e- 


17 


9.0e- 


■16± 2.9e-17 


3.4± 0.1 


4.5± 0.2 


4.0± 0.4 


4.0+ 0.1 


35 




9.7e- 


■16± 6.7e- 


■17 


4.8e- 


■16± 3.5e- 


17 


1.7e-16± 2.9e- 


17 


1.2e- 


■15± 3.3e-17 


5.1±0.2 


5.1±0.2 


4.9± 0.4 


4.9+ 0.1 


36 




6.7e- 


■16± 4.4e- 


■17 


2.2e- 


■16± 2.2e- 


17 


9.9e-17± 3.0e- 


17 


9.3e- 


■16± 2.9e-17 


4.0± 0.1 


2.3±0.1 


3.3± 0.5 


4.4+ 0.1 



Table 3. Measurements described in Sec. 13.21 for the apertures defined on the Bry image and the H2 image. The Une fluxes are in 
erg s"' cm"^, and the line Gaussian <j are in A. 



aperture 


name 


2-lS(2) flux 


l-OS(O) flux 


2-lS(l)flux 




0-2-18(2) 


a l-OS(O) 


o-2-lS(l) 


3 




2.0e-16+ 3.5e- 


■17 


5.7e-16+ 3.7e- 


17 


4.0e-16+ 4.9e- 


17 


5.2+ 0.5 


4.9+ 0.2 


5.2+ 0.3 


4 




2.0e-16+ 2.9e- 


17 


4.8e-16+ 3.0e- 


17 


3.3e-16+ 5.8e- 


17 


3.8+ 0.3 


3.5+0.1 


4.7+ 0.4 


7 




< 2.2e-16 




3.4e-16+ 2.3e- 


17 


2.0e-16+ 3.2e- 


17 




3.3+0.1 


2.8+ 0.2 


8 




1.6e-16+ 4.5e- 


17 


7.3e-16+ 3.6e- 


17 


4.3e-16+ 4.8e- 


17 


5.1+ 0.7 


4.3+0.1 


4.3+ 0.2 


9 


M4 


1.5e-16+ 3.8e- 


17 


7.0e-16+ 4.2e- 


17 


3.9e-16+ 5.4e- 


17 


4.5+ 0.6 


4.6+ 0.1 


4.1+0.3 


11 




2.1e-16+ 4.6e- 


17 


3.5e-16+ 3.4e- 


17 


2.0e-16+ 2.4e- 


17 


7.1+0.8 


4.1+ 0.2 


2.5+0.1 


13 




< 3.4e-16 




< 1.6e-16 




< 2.7e-16 










15 




1.9e-16+ 2.6e- 


17 


6.6e-16+ 3.9e- 


17 


4.5e-16+ 3.6e- 


17 


2.8+ 0.2 


4.2+ 0.1 


4.4+ 0.2 


16 


M5-4 


7.0e-17+ 2.2e- 


17 


3.9e-16+3.6e- 


17 


2.2e-16+ 5.1e- 


17 


2.0+ 0.3 


4.3+ 0.2 


5.0+ 0.6 


18 


M5 


< 1.9e-14 




5.5e-16+ 4.6e- 


17 


4.0e-16+ 6.0e- 


17 




3.8+ 0.2 


3.3+0.3 


21 


M5-2 


< 2.2e-14 




5.7e-16+ 5.2e- 


17 


2.5e-16+ 4.1e- 


17 




4.9+ 0.2 


2.8+ 0.2 


22 


M5-3 


1.3e-16+ 3.2e- 


17 


6.7e-16+ 4.0e- 


17 


3.6e-16+ 3.8e- 


17 


3.4+ 0.4 


4.3+ 0.1 


3.6+ 0.2 


23 




1.6e-16+ 2.9e- 


17 


6.1e-16+ 3.7e- 


17 


3.5e-16+ 3.7e- 


17 


3.2+ 0.3 


4.0+ 0.1 


3.5+ 0.2 


28 




6.7e-17+ 1.9e- 


17 


2.2e-16+ 2.9e- 


17 


1.8e-16+ 3.2e- 


17 


1.8+ 0.3 


3.6+ 0.2 


2.9+ 0.3 


29 




< 3.6e-14 




5.7e-16+ 7.3e- 


17 


3.0e-16+ 5.1e- 


17 




4.8+ 0.3 


3.9+ 0.3 


31 


M6 


< 3.7e-14 




5.1e-16+ 8.1e- 


17 


3.2e-16+ 5.1e- 


17 




4.4+ 0.4 


3.6+ 0.3 


33 




< 2.2e-14 




5.2e-16+ 4.2e- 


17 


2.8e-16+ 4.8e- 


17 




5.9+ 0.2 


5.7+ 0.5 


34 




1.2e-16+ 2.9e- 


17 


3.3e-16+ 2.8e- 


17 


3.3e-16+ 5.4e- 


17 


4.6+ 0.6 


3.3+ 0.1 


6.8+ 0.6 


35 




< l.le-16 




4.2e-16+ 3.4e- 


17 


1.2e-16+ 2.0e- 


17 




6.4+ 0.3 


1.7+ 0.1 


36 




1.4e-16+ 2.9e- 


17 


3.0e-16+ 2.8e- 


17 


1.4e-16+ 2.0e- 


17 


5.2+ 0.5 


3.4+ 0.2 


2.4+ 0.2 



Table 4. Measurements described in Sec. 13.21 for the apertures defined on the Bry image and the H2 image. The line fluxes are in 
erg s"' cm"^, and the line Gaussian cr are in A. 
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